Internal loss of superconducting resonators induced by interacting two level systems 
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In a number of recent experiments with microwave high quality superconducting coplanar waveg- 
uide (CPW) resonators an anomalously weak power dependence of the quality factor has been 
observed. We argue that this observation implies that the monochromatic radiation does not sat- 
urate the Two Level Systems (TLS) located at the interface oxide surfaces of the resonator and 
suggests the importance of their interactions. We estimate the microwave loss due to interacting 
TLS and show that the interactions between TLS lead to a drift of their energies that result in a 
much slower, logarithmic dependence of their absorption on the radiation power in agreement with 
the data. 



High quality superconducting CPW resonators are 
used in a number of diverse fields, ranging from as- 
tronomical photon detection [l], 0] to circuit quantum 
electrodynamics In these applications, the CPW 

resonator is operated in a regime of low temperature 
(~ lOmk) and low excitation power (single photon). The 
performance of these devices is directly related to the 
resonator quality factor, Q, defined by the photon de- 
cay rate, j p h, as Q = u>o/"/ p hi where cjq is the resonance 
frequency. j p h is given by the sum of the escape rate 
from the resonator and the rate of the intrinsic decay; 
the latter sets the limit on the performance. At low tem- 
perature and single photon regime, the intrinsic decay is 
usually attributed to the excitations of TLS located in 
the dielectrics (bulk and surfaces) surrounding the res- 
onator Q. This belief is strongly supported by the ob- 
servation of temperature-dependent resonance frequency 
shift that closely agrees with the one predicted by the 
conventional theory of microwave absorption of TLS in 
glasses [H 0] ■ According to this theory, one expects also 
that, as the power of the radiation applied to resonator is 
increased, TLS in the dielectrics get saturated, thereby 
limiting the maximal power that can be dissipated by 
photons. This results in a strong power dependence of 
the quality factor: Q cx \[P above a critical power P c . 
This power dependence is indeed observed in many res- 
onators characterized by intrinsic loss tangent ~ 10~ 3 at 
very low powers [ToL ITll ] . However resonators character- 
ized by lower intrinsic loss at low powers typically show 
much weaker power dependence 12- lfjj|. 

In fact, careful fitting of the loss versus power to the 
empirical equation Q cx (1 + P/P c ) v gives <p ~ 0.03—0.16 
for the resonators made of Nb and Al on sapphire and 
of Nb on undoped silicon (l3j ; similarly high-Q single- 
layer resonators of different geometries made of Al on 
sapphire substrate show <p ~ 0.1 — 0.2 [15j. Corrections 
of the conventional theory of TLS which include the ge- 
ometric dependence of the applied electric field failed to 
reproduce the weak power dependence observed experi- 



mentally [I2|,llj-[16|. 

The failure of the conventional theory of TLS to predict 
the power dependence of the quality factor for the high 
quality resonators is an indication of a serious gap in our 
understanding of TLS in amorphous insulators. In this 
Letter we argue that in high-Q superconducting CPW 
resonators the TLS located in the interface oxide surfaces 
are subject to stronger interactions than the TLS located 
in the bulk dielectrics. These TLS interactions lead to a 
drift of the TLS energies that results in a logarithmic 
dependence of their absorption on the radiation power, 
P, in agreement with the data. We begin by sketching 
the most important assumptions that lead to the Q cx 
\[P prediction of the conventional theory of microwave 
absorption of TLS; then we discuss the implications of the 
much weaker power dependence reported experimentally 
in high-Q superconducting CPW resonators. 

In the conventional theory, TLS are described by 
pseudo-spin operators, S, and are characterized by the 
uniform distribution of the energy difference, e, between 
their ground and excited state. In the basis of the eigen- 
states the Hamiltonian has a simple form H = eS z . The 
ground and the first excited state of the TLS correspond 
to the quantum superposition of the states characterized 
by different atomic configurations. Each TLS is charac- 
terized by a dipole moment p = p (sin 6S X + cos9S z ), 
which is an operator with both diagonal and off-diagonal 
components, p denotes the difference between the dipole 
moments in the two different atomic configurations, its 
magnitude po = |p| sets the scale of the dipole mo- 
ment. describes the fact that the eigenstates of the 
dipole correspond to the superposition of its states in 
real space. Because many dipoles have exponentially 
small amplitude to tunnel between different positions in 
real space, the parameters and e are assumed to have 
distribution V(e,0)dedO ~ v/OdedO for small 6, where 
v = 10 20 /cm 3 eV is the typical density of states of TLS 
1 71 ] . The interaction between different TLS is essentially 
of a dipolc-dipolc nature with an effective strength given 



2 



by the dimensionless parameter A = fpg/e, here e is the 
dielectric constant of the medium, it can also be viewed as 
coming from TLS coupling to elastic strain 18] . Straight- 
forward analysis shows that the same parameter also con- 
trols the phonon mean free path at low temperatures [l9| . 
The direct measurement gives values of A f=s 1CP 3 in bulk 
materials, so the interaction between TLS is usually as- 
sumed to be irrelevant. 

The intrinsic microwave loss is due to the coupling of 
the electric dipole moment p of the TLS to the applied 
electric field Ecoswt of the resonator. The resonator qual- 
ity factor is related to the imaginary part of the dielectric 
function, e(uj), i.e. oc |3e|. At microwave frequen- 

cies and low temperatures, only the resonant contribution 
to the electric susceptibility tensor is relevant. One can 
compute the change in the dielectric function due to the 
resonant response of an ensemble of TLS by assuming 
that TLS are relaxed by phonons and using the Bloch 
equations that neglect the interaction between TLS [3] ■ 
In fact, one can get the power dependence of quality fac- 
tor Q by following a more general, qualitative argument. 
In the steady state driven by sinusoidal electric field, the 
power dissipation density is given by the time average 
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where D= cE is the displacement vector. In the absence 
of TLS-TLS interactions, each TLS can decay from its 
excited state only by emitting a phonon; the rate of this 
emission gives the decay rate, L(e), of its excited state. 
The maximal power that a given TLS can absorb from 
the photons in the resonator is sT. At small applied pow- 
ers, the TLS that absorb photons have energies close to 
the frequency, wo, of the photon field, i.e. |e — cjo| < L. 
It is convenient to characterize the effect of the electric 
field on TLS by their Rabi frequency ft = ip • E sin 0; for 
small ft < T TLS at resonance are excited with proba- 
bility (f2/L) 2 . Since their density is vT, the total power 
dissipation density reads (Pd) = ve^l 2 /T. Different is the 
situation at larger applied powers. There, even the TLS 
with energies further away from the frequency of the res- 
onator, wo, get excited. For TLS exactly at resonance, 
the electric field E cos oj^t results in oscillations with fre- 
quency il. When f2 > T, TLS with energies \e — ujq\ <il 
are excited. Since their density is z^fi, the total dissipated 
power density reads (Pd) — veTfl. By substituting it into 
Eq. [TJ one recovers immediately the quality factor power 
dependence Q cx VP. 

The generality of these qualitative arguments implies 
that, in the resonators where a much weaker power de- 
pendence of Q has been observed, the dipoles that absorb 
the radiation do not get saturated as the applied power P 
is increased. This is possible if the interactions between 
dipoles are not negligible and lead to energy diffusion. 

The bulk of evidence indicates that the loss in high- 
Q superconducting CPW resonators is due to TLS lo- 
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Figure 1: Drift of the energy of a given TLS in the field of a 
few fiuctuators. 



cated on the interface oxide surfaces (metal-air, metal- 
substrate or substrate-air) even though it 
remains unclear which of the surfaces is more relevant{20] . 



The estimates show [21| that the concentration of TLS 
in these thin surfaces is higher than in the bulk and con- 
sequently the average interaction between TLS is larger, 
implying that one has to develop the theory of microwave 
absorption in interacting TLS in these resonators. In 
contrast, the resonators that show Q oc VP" and low in- 
trinsic loss in single photon regime are made of Nb on 
Si0 2 /Si [HI and AlO^ coated Nb on Si substrate [n| 
and thus contain a significant amount of bulk amorphous 
dielectrics (SiC>2 or AlO^). It is natural to assume that 
the large intrinsic loss in these resonators is due to TLS 
located in the dielectric bulk which is described well by 
the conventional theory of independent TLS. 

Developing the full theory of interacting TLS is a very 
difficult problem that was first discussed by Yu et al. 22 1 
and still remains controversial (23l426| . Fortunately, as we 
shall see below, one does not need to solve the full theory 
of interacting TLS in order to estimate the internal loss 
of superconducting resonators due to interacting TLS. 
In our approach, we assume that the effective degrees 
of freedom that remain active in the amorphous oxides 
surrounding the superconductors are described by fluc- 
tuating dipoles: some of these dipoles are characterized 
by fast transitions between their states (esin^ ~ a;) and 
relatively small decoherence rates, i.e. r <C ui, and are 
effectively coherent TLS; other dipoles are instead slow 
and characterized by decoherence times shorter than the 
typical time between the transitions, we shall call them 
fiuctuators. Due to the interaction between fiuctuators 
and TLS, the frequency of the TLS jumps when the fiuc- 
tuators in its vicinity change its state. This translates 
into the fact that instead of staying constant the energy 
of a given TLS drifts with time (see Fig. [IJ. Notice that 
this phenomenological model is in a full agreement with 
experiments on char ge n oise performed in superconduct- 
ing SET and qubits [27j]. 

We now use this phenomenological model to estimate 
the internal loss of the CPW resonators. We begin with 
qualitative arguments. Due to the interaction with fiuc- 
tuators, the TLS with energy level e in resonance with the 
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applied electric field stays effectively in resonance only for 
a short dwell time r = j^ 1 , where 7 is the combined re- 
laxation rate of the fluctuators that affect it (see FigfTJ). 
If the time r is short, i.e. rft <C 1, the TLS gets into 
the excited state with probability V c ~ (Qi~) 2 and then 
dissipates the energy away from the resonance. The av- 
erage dissipation rate of this process is f = jV e (x £l 2 /"i- 
Only TLS with energy level located within energy 7 
from the resonator frequency u> contribute to this pro- 
cess: their density is and the total dissipated power 
density reads (Pd) = veil 2 JJ max "P(7)e?7 , where V(j) = 
Po /"f is the probability distribution of the fluctuator re- 
laxation rates [28]. After integration, we obtain that 
(Pd) = vPqeSI 2 In ("fmax/Q) and by inserting it into EqQ] 
we find that in this case Q oc \n(il/-f max ) oc InP, i.e. 
there is only a slow, logarithmic dependence of the qual- 
ity factor with the applied power P. 

These qualitative arguments can be confirmed by more 
quantitative analytical computations. Before proceeding 
with the analytic derivation, we stress that the essen- 
tial ingredient of the phenomenological model is the large 
value of the jumps in the TLS frequency, Set ^ T caused 
by fluctuators. This implies that the interaction between 

the high frequency dipole and the fluctuators in its vicin- 

2 

ity is large: V(r) = >■ T. Because the typical dis- 
tance between thermally activated TLS is f = p -1 / 3 , with 
p = vT, the typical interaction is V(f) — XT. Compar- 
ing it with the relaxation rates of T ~ 10 6 s _1 observed 
experimentally at frequencies uj ~ 10 GHz (and expected 
theoretically for phonon emissions) we see that in typi- 
cal experimental conditions de ~ V for A ~ 10~ 3 . The 
condition Se 3> T, however, can be satisfied for lower fre- 
quencies or in materials with larger A. The estimates of 
the dissipation per unit volume of surface oxides show 
that there their density of TLS is at least 10 times larger 
than in the bulk amorphous materials [2lj. This trans- 
lates into 10 times larger value of the parameter A and 
a completely different physics of microwave dissipation 
described above. 

We now give the details for the analytical derivation 
of the Q oc y/P and its generalization for the model that 
takes into account energy drifts of TLS caused by fluc- 
tuators. The collection of TLS is characterized by the 
Hamiltonian H int = £\ e { Sf + \ J2 itj V ab (ri - r 3 )p a p b . 
Here V a b{r) is the interaction between TLS that is due to 
their dipole moments or virtual phonon exchange. In ei- 
ther case, in the static limit, the interaction scales as 
V(r) r~j 1/r 3 . At zero temperatures and for a small 
value of the parameter A <C 1, the Hamiltonian Hi nt 
gives a coherent dynamics of the individual TLS. In this 
limit, the levels at high frequencies are broadened by 
phonon emission [i| that leads to the relaxation rate 
Fi oc e 3 . At non-zero temperatures the levels are fur- 
ther broadened by dephasing caused by thermal phonons 
and other TLS [H, So|. The combined effect of all 



these processes on a single TLS can be described by 
the Bloch equations: f t (S(t)) = [(S(t)) x B(t)} - R(t), 
where R(t) = (T 2 {S x ),T2{Sy),Ti({S g ) - m)) is the re- 
laxation matrix, B(£) is the total field acting on the 

TLS and m = — tanh [e/2T] denotes the equilibrium pop- 
ulation of the TLS levels and T is the temperature. 
In the conventional theory of TLS, this field has two 
components: B = B° + B 1 , with B° = (0,0, e) and 
B 1 = (O,0,ft')coswi , where ft' = §pEcos0. The pres- 
ence of the fluctuators results in the additional time de- 
pendence of the energy levels, e t = e + Here 
denotes a multi-level telegraph noise signal characterized 
by switching rate 7. 

The nonlinear Bloch equations are dramatically sim- 
plified in realistic conditions, because the feasible elec- 
tric field acting on the TLS has very little effect on 
them away from the resonance. Formally, in the sta- 
tionary solution of the Bloch equations driven with 
Ecoswt, all spin components oscillate with frequencies 
nu)\ S = ^2 S n exp(— inuit) with S„ — S*L n . Relax- 
ation to the stationary solution is determined by the de- 
cay rate ri -C T2 <C oj which sets the longest time scale 
in the problem. It can be described by allowing slow time 
dependence of S„ components. Finally, introducing the 
operators 5* = S x ± iS v and leaving only the leading 
resonant terms the Bloch equations reduce to: 

= nS z Q - (w - s t + iT 2 )S+ (2) 

^ = ^ 1 + -r 1 (5 2 °-m) (3) 

The macroscopic response is given by the average po- 
larization, P w of TLS at frequency to: 

P u = I (p sin 0St) (4) 

where the average is taken over the distribution of TLS 
and their dipole moments. The quality factor is directly 
related to P w : Q" 1 oc |9fP w | /E. 

If the jumps induced by the fluctuators are small, 
i.e. £ < T2, the stationary solution of Eqs. (|2l3[) gives the 
spin component: 

+ fl(uj-e- iT 2 )m 

1 (w- e )2 + rf + 02r 2 r 1 - 1 ' w 

At small fields, one can neglect the last term in the de- 
nominator of Eq.([5]). By substituting Eq.© into Eq.((3]) 
and then averaging over the distribution of TLS, one gets: 

|3PJ = ^(sin 2 0)AE (6) 
6 

and consequently a Q factor that does not depend on the 
field strength E. 

At large fields (ft 2 > ri^), the third term in denomi- 
nator of Eq.© dominates the second and this results in 
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a rapid decrease of the response. By averaging over the 
distribution of TLS, one gets: 



TTin 
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which translates into the usual Q ~ V P dependence. 

In the opposite case of large jumps £ 3> L2 one should 
solve the full dynamical Eqs. (|2|3p with noise £(t). The 
details of the solution depend on the relation between the 
relaxation rates Fi,F2 and the jump rate 7. We discuss 
first the simplest case of large 7 in which a given TLS 
spends most of its time away from the resonance, so that 
when it moves back into the resonance its magnetization 
Sq and S± have their equilibrium values: Sq = m and 
= 0. In this case the only quantity that controls the 
response is the rate, 7, of incoming and outgoing jumps 
into the resonance. The stochastic nature of this process 
implies that, in order to find the average value of 
that determines the response to electric field, we have 
to solve the equation for the probability, g(s), to find 
the resonant TLS characterized by the three dimensional 
vector s = (x, y, z), where x and y are the real and the 
imaginary parts of S^ and z = Sq. This probability 
obeys the evolution equation: 



do 
di 



d ( ds 



ds \ dt 



g ) = 7 [8{z - m)5(x)8(y) - g] 



where ds/dt is given by Eqs. ([2l3p with constant e. The 
equation is further simplified in the limit of large 7 > 
fl 3> Li, L2 where the analytic solution gives (see Supple- 
mentary material), after averaging over the distributions 
of e and 7, "P(e,7) = vPq/^: 
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which results in a weak, logarithmic dependence of the 
quality factor Q on the strength of the electric field. 

A similar logarithmic dependence occurs in the cases 
when the dephasing rate is larger than the jump rates 
as well. In this case, one can neglect the time derivative 
terms in Eq. (jSJ and express S± through Sq and get the 
general evolution equation: 

%^ + TT [(ri(z - m) - E k z) g k ] = JknQn (9) 
at oz 

where gk(z,t) is the probability for a given TLS to have 
value Sq = z while subjected to the effective driving field 
Sfe = 2 r 2 /[(o; — £fc) 2 + rf] in the state k of the fluctua- 
tors. The matrix jki is the matrix of the transition rates 
between the states of the fluctuators. Similarly to the 
case of the small dephasing rate the saturation of TLS 
does not happen if they stay mostly away from the res- 
onance, so that the probability to find a given TLS in 
resonance n < T±/j (see Supplementary material). Be- 
cause the number of states of classical fluctuators grows 



exponentially with their number this condition is satisfied 
even for a moderate number of fluctuators that affect a 
given TLS. In this situation the dissipation is dominated 



(7) by TLS with 7 > 5 - 2 /F 2 



I3PJ = 
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(10) 



Both Eqs ()8|10|) lead to the weak, logarithmic depen- 
dence of Q at large electric fields (J7 2 > L2L1), while at 
smaller fields, the imaginary part of the average polar- 
ization is given by Eq.([B]) and Q is constant. 

As a final remark let us notice that in this model the di- 
mensionless coupling between TLS remains small, A <C 1. 
This implies that thermodynamics properties and the real 
part of the dielectric constant are not affected by the in- 
teraction. Therefore we expect the same temperature- 
dependent resonance frequency shift as predicted by the 
conventional theory of TLS in agreement with experi- 
ments |8j, l9(. 

In conclusion, we have shown that the resonance ab- 
sorption from TLS is strongly affected by their interac- 
tion with classical fluctuators; this does not change the 
absorption at low powers but changes the square root 
dependence of the absorption into the logarithmic one 
at high powers. The important condition is that each 
fluctuator should move the TLS frequency more than its 
width due to the relaxation rate. This translates into 
a higher (~ 10 times larger) concentration of TLS than 
the typical density in amorphous bulk materials.. We in- 
terpret the results of recent experiments displaying slow 
power dependence of the quality factor in high-Q super- 
conducting CPW resonators as the evidence for a large 
concentration of TLS located in the interface oxide sur- 
faces of the resonator. Very likely it implies that these 
TLS have a different nature, e.g. localized electron states 
at the superconductor oxide boundary. 

Note added: recent experiment directly observed the 
energy drift of TLS which is the basis of our model 31 1 
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